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Abstract 

The addition of non-renormalizable terms involving the Higgs fields to the MSSM (BMSSM) 
ameliorates the little hierarchy problem of the MSSM. For neutralino dark matter, new regions 
for which the relic abundance of the LSP is consistent with WMAP (as the bulk region and the 
stop coannihilation region) are now permitted. In this framework, we analyze in detail the direct 
dark matter detection prospects in a XENON- like experiment. On the other hand, we study 
the capability of detecting gamma-rays, antiprotons and positrons produced in the annihilation 
of neutralino LSPs in the Fermi and oncoming AMS-02 experiments. 
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1 Introduction 



The smallness of the quartic Higgs couphng in the framework of the minimal supersymmetric 
standard model (MSSM) poses a problem. The tree level bound on the Higgs mass is violated, 
and large enough loop corrections to satisfy the lower bound on the Higgs mass suggest that the 
stop sector has rather peculiar features: at least one of the stop mass eigenstates should be rather 
heavy and/or left-right-stop mixing should be substantial [1]. 

The situation is different if the quartic Higgs couplings are affected by new physics. If the new 
physics appears at an energy scale that is somewhat higher than the electroweak breaking scale, 
then its effects can be parametrized by non-renormalizable terms. The leading non-renormalizable 
terms that modify the quartic couplings are [2l[3llll[5l|6l[71[8l|9l [TU] : 

VFbmssm = ^{HuHdf + ^Z{HuHdf, (1) 
where 2 is a SUSY-breaking spurion: 

Z = e'^nisusy (2) 

The first term in equation ([1]) is supersymmetric, while the second breaks supersymmetry (SUSY). 
In the scalar potential, the following quartic terms are generated: 

2eiHuHd{HlHu + H^Ha) + e2{HuHdf, (3) 

where 

" IvT' Ai- 

However, let us note that these operators Beyond the MSSM (BMSSM) may destabilize the 
scalar potential. If 4|ei| > €2, the effective quartic coupling along one of the D-flat directions is 
negative, causing a remote vacuum to form in the presence of which the electroweak vacuum could 
become metastable 

The interplay between the Higgs sector, the stop sector, and the non-renormalizable (NR) ope- 
rators has interesting consequences for the MSSM baryogenesis |12| I13| . The window for MSSM 
baryogenesis is extended and, more importantly, can be made significantly more natural. In addi- 
tion, these operators have implications for yet another cosmological issue, and that is Dark Matter 
(DM) |14| [T5l 113) . The phenomenology of the BMSSM at colliders has been studied in reference 
|16] : implications for fine-tuning have been analyzed in reference [17j. 

One of the attractive features of the MSSM is the fact that the lightest R-parity-odd parti- 
cle (LSP) is a natural candidate for being the dark matter particle. Progress in experimentally 
constraining the MSSM parameter space restricts, however, the regions where the dark matter is 
quantitatively accounted for to rather special regions of the MSSM: the focus point region, with 
surprisingly heavy sfermions; the funnel region, where the mass of the CP-odd neutral Higgs scalar 
is very close to twice the mass of the LSP; the co-annihilation region, where the mass of the scalar 
partner of the right-handed tau is very close to the mass of the LSP; and the bulk region, where 
the bino-LSP and the sleptons are light. 
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The effects of the NR operators are potentially important for two of these four regions. First, 
these operators give rise to a new Higgs-Higgs-higgsino-higgsino interaction Lagrangian, 



/X* 



2{HuHd){HuHd) + 2{HuHd){HuHd) + {HMiHuH^) + (HMiHuHd) + h.c, (5) 



which contributes to the annihilation process of two higgsinos to two Higgs particles. This effect is 
relevant when the dark matter particle has a significant component of higgsinos, as is the case in 
the focus point region. Second, as mentioned above, these operators modify the relation between 
the light Higgs mass and the stop masses. This effect can be important in the bulk region within 
models where the slepton and stop masses are related, such as the mSUGRA models. 

The effect of these operators on the relic density was studied in detail in references |15| [T3] . It 
was found that new regions yielding the correct relic density can arise, especially in the bulk and 
co-annihilation region. In this work, we examine how the detection prospects of the MSSM are 
modified by the introduction of terms in equation ([3]) in the superpotential, focusing ourselves on 
the re-opened regions. 

We shall evaluate the detection perspectives for the main four kinds of detection usually consid- 
ered in the literature, namely direct detection in a XENON-like experiment, gamma-ray detection 
from dark matter annihilations in the galactic center for the Fermi mission |18) as well as positron 
and antiproton detection coming again from dark matter annihilations in the galactic halo and for 
the oncoming AMS-02 |19) experiment. 



2 The model 

The BMSSM framework, if relevant to the little hierarchy problem that arises from the lower bound 
on the Higgs mass, assumes a new physics scale at a few TeV. Since the new degrees of freedom 
at this scale are not specified, the effect of the new threshold on the running of parameters from a 
much higher scale cannot be rigorously taken into account. It therefore only makes sense to study 
the BMSSM effects in a framework specified at low energy. In order to demonstrate some of the 
most interesting consequences of the BMSSM operators for dark matter, we shall employ the two 
sets of parameters explored in reference |13) : a model where all sfermion masses are correlated, 
and a model where the only light sfermions are the stops. The first model demonstrates how the 
so-called bulk region is re-opened, even for correlated stop and slepton masses. The second model 
incorporates the interesting process of stop co-annihilation. For both models we focus our attention 
mainly on regions where the stops are light, since the main motivation for the BMSSM operators is 
to avoid a heavy stop (which is the cause of the little hierarchy problem). Previous analysis in the 
context of the MSSM with a light stop was done in references |20| I21j . 

Within this framework, we calculate the dark matter relic density, and the direct and indirect 
detection prospects in the presence of the new ej couplings. We used a modified version of the code 
micrOMEGAs |22| I23j . where we implemented the BMSSM Higgs-Higgs-higgsino-higgsino couplings 
of equation in order to calculate the relic density as well as the cross-sections and decay 
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channels relevant for dark matter detection. The leading ej-induced corrections to the spectrum, 
were implemented using the code SuSpect ["i^. 

2.1 Correlated stop-slepton masses 

The first scenario considered contains correlated stop and slepton masses, just as the most studied 
MSSM scenarios, such as the mSUGRA [251 [261 [271 [28] or cMSSM frameworks. In this case, the 
neutralino LSP is a bino-like state annihilating to the standard model leptons via light slepton 
exchange. However, this scenario, known as the 'bulk region', is highly constrained due to the 
experimental lower bound on the Higgs mass. In general, in order to fulfill such a constraint either 
heavy or mixed stops are required [29] . In the BMSSM, nevertheless, it is possible to re-open the 
bulk region regardless of the structure of the stop sector. 

In order to allow for a simple comparison with mSUGRA- models, we focus the attention to the 
parameters 

tan/3, mi/2, mo, Aq, sign(/u). (6) 

Let us emphasize again that one should not think about this set of parameters as coming from 
an extended mSUGRA model, since the effects of the BMSSM physics at the few TeV scale on 
the running are not (and cannot) be taken into account. In addition, we have two extra BMSSM 
parameters: ei and €2- 

In practice, we make discrete choices of tan/3, Aq, sign(/i), ei and €2, and scan over ttt-q and 
mi/2- We focus our attention on moderate values of mi/2 and mo because we are mainly interested 
in light sfermions and the bulk region. We also use Aq = GeV and /i > in the whole analysis. 

As pointed out in reference |13j. we would like to emphasize several points regarding the present 
scenario: 

1. A generic point in the former parameter space usually gives rise to a too low annihilation 
cross-section of neutralino LSP and hence to a too large relic density, in confiict with the 
WMAP measurements |30j. 

2. However, for moderate mi/2 and low tuq values, there is a region where the LSP is almost 
degenerate in mass with the lightest stau (fi), enhancing the co-annihilation cross-section 
X?-r. 

3. Another region giving rise to relic density in agreement with the WMAP measurements ap- 
pears for mi/2 ~ 120 GeV. This is the '/i-pole' and the 'Z-pole' region in which nih ~ mz ~ 
2m^o, and the s-channel Higgs and Z boson exchange is nearly resonant, allowing the neu- 
tralinos to annihilate efficiently |31) . Let us note that for the ordinary mSUGRA case, this 
region is already excluded by LEP measurements. 

4. For negative enough ei values, the uplift of the Higgs mass generates a splitting among the 
'/i-pole' and the 'Z-pole' regions, with the former now evading LEP constraints. 
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The latter point concerning the Higgs boson mass is the most significant effect of the BMSSM 
operators. Within the MSSM with mSUGRA-hke correlations, the bound on the Higgs mass strongly 
constrains mi/2- In contrast, in the presence of ei = O{—0.1), the full region for which the correct 
value of the relic abundance is obtained is allowed. Let us emphasize that in the plain MSSM 
scenario at stake the bulk region is already ruled out because it gives rise to a too light Higgs boson, 
in contradiction to LEP2 data |32j . However, the introduction of the NR operators in equation ([1]) 
uplifts the latter mass up to rufi > 130 GeV or rufi > 155 GeV for tan /? = 10 or 3 respectively. 

In the mo region that we are considering here, the impact of the BMSSM operators on the mass 
of the neutralino LSP is rather limited. The reason is that in the bulk region the LSP is mostly 
bino-like, while the BMSSM operators affect the higgsino parameters. 

Concerning precision electroweak data and low energy processes, it is important to realize that 
the new physics that generates the non-renormalizable operators can directly modify the constraints 
that come from these measurements. Ignoring this point, it is still possible to identify regions in 
the parameter space favored by the WMAP data which satisfy all such low energy constraints. The 
relevance of the BMSSM lies in the fact that constraints involving the Higgs are decoupled from 
constraints involving the stop sector. 

2.2 Light stops, heavy sleptons 

In order to continue with the analysis of scenarios with light unmixed stops, we focus on a set of 
low energy parameters very different from the previous subsection. Explicitly, in addition to the 
BMSSM €i parameters, we consider the following set of parameters: 

Ml, /X, tan/3, Xt, mu, mq, mj, tua, (7) 

where mj is a common mass for the sleptons, the first and second generation squarks, and We 
further use Mi = | tan^ 9w M2 ~ ^ M2. To demonstrate our main points, we fix the values of all 
but two parameters as follows: ei = or —0.1, €2 = 01 +0.05, tan/3 = 3 or 10, Xt = 0, mu = 210 
GeV, niQ = 400 GeV, = = 500 GeV. This scenario gives rise to relatively light stops: 

< 150 GeV, 370 GeV < m^^ < 400 GeV. (8) 

We scan over the remaining two parameters. Mi and fi. 

Again, as pointed out in reference |13) . in the prescribed framework one can identify four regions 
in which the WMAP constraint is fulfilled: 

1. The 'Z-pole' region in which the LSP is very light, m^o ~ ^Mz ~ 45 GeV, and the s-channel 
Z exchange is nearly resonant. This region is not ruled out only in scenarios where the mass 
splitting between Mi and M2 at the electroweak scale is very large. 

2. The '/i-pole' region in which the LSP is rather light, m^o ~ ^M/j, and the s-channel h exchange 
is nearly resonant, allowing the neutralinos to annihilate efficiently |31) . 
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3. The 'mixed region' in which the LSP is a higgsino-bino mixture [33], Mi ~ fi, which enhances 
(but not too much) its annihilation cross-sections into final states containing gauge and/or 
Higgs bosons: XiXi ^ W+W , ZZ, Zh and hh. 

4. The 'stop co-annihilation' region, in which the LSP is almost degenerate in mass with the 
lightest stop (ti). Such a scenario leads to an enhanced annihilation of sparticles since the 
Xi — ii co-annihilation cross-section [341 [35l [36] is much larger than that of the LSP. 

Let us first consider the case where ei = e2 = 0. The region at Mi ~ mzf^ ~ fnh/'^ corresponds 
to the s-channel exchange of an almost on-shell Higgs or Z boson. Note that when 2 m^o is too close 
to the Higgs or Z mass pole, the LSP annihilation is too efficient and leads to a much too small 
^DM In any case, for the Higgs mass values obtained here, nih ~ 85(98) GeV for tan/3 = 3(10), 
this region is already excluded by the negative searches for chargino pairs at LEP2 [32] . 

The region close to M2 ^ ^ ^ 200 GeV corresponds to the LSP being a bino-higgsino mixture 
with sizeable couplings toW,Z and Higgs bosons, allowing for reasonably large rates for neutralino 
annihilation into XiXi ~^ W~^W^, ZZ, hZ and hh final states. Above and below the band, the 
LSP couplings to the various final states are either too strong or too weak to generate the relevant 
relic density. 

Finally, for larger fi values, the mass of the lightest neutralino approaches the mass of the 
lightest stop leading to an enhanced co-annihilation cross- section: Xi^i W+b, gt 90%). 
Also, to a lesser extent (~ 5%), the annihilation cross-section of the stop NLSP contributes to the 
total cross-section by the process ti ti — )• 5 g. 

Next we consider the ei = —0.1 case. The features of the DM allowed regions are similar to the 
previous case. The main difference comes from the important enhancement of the Higgs mass due 
to the presence of the BMSSM operators. In this case it is possible to disentangle the Z and the 
h peaks, since the Higgs- related peak moves to higher M2 values, due to the increase of the Higgs 
mass: nih = 122(150) GeV for tan/3 = 10(3). Furthermore, the latter peak is no longer excluded 
by chargino searches. 



3 'Detectability' definition 

Let us now define what we shall be meaning by saying that a certain parameter space point is 
detectable. We employ a method based on the quantity. Consider whichever mode of detection: 
direct or indirect in any of the three channels (7-rays, e^,p) we shall be considering. In all three 
modes, what is finally measured is a number of events per energy bin. Let us call A^**^ the number 
of signal (dark matter - induced) events in the i-th. bin, the nature of which depends on the specific 
experiment, N^^^ the corresponding background events in the same bin, and Nf°^ the sum of the 
two. The variance xj in every bin is defined as: 

^tot ■ v-"; 
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Then, the condition that we impose to characterize a point as detectable, is that at least in one 
energy bin xf ^ 5.8. In Gaussian error terms, this corresponds to a 95% CL. 

4 Direct detection 
4.1 Differential event rate 

In spite of the experimental challenges, a number of efforts worldwide are actively pursuing to 
directly detect WIMPs with a variety of targets and approaches (for previous works see, e.g. refe- 
rences \37 \ 1381 [39] ). Many direct dark matter detection experiments are now either operating or in 
preparation. All these experiments measure the number of elastic collisions between WIMPs and 
target nuclei in a detector, per unit detector mass and per unit of time, as a function of the nuclear 
recoil energy Er- The detection rate in a detector depends on the density po — 0.385 GeV cm~^ 
|40j and velocity distribution /(%) of WIMPs near the Earth. We assume a Maxwellian halo for 
WIMP's velocity in the rest frame of our galaxy (for a recent treatment of non-Maxwellian case, see 
e.g. reference |41|). taking into account the orbital motion of the solar system around the galaxy, 
and neglecting the motion of the Earth around the Sun |42| : 

ffy ) = -^^^ L-K-1.05.o)V-§ _ g-K+1.05.o)V^'gl (10) 

where vq ~ 220 km/s is the orbital speed of the Sun around the galactic center. In general, the 
differential event rate per unit detector mass and per unit of time can be written as: 

dN aopo r fi^x 

dEf 2 m,^ m 



F{ErY -^dv^, (11) 



where ctq is related to the WIMP-nucleon cross-section, cr^-p, by do = cr^-p ■ (Anir/Mr)'^ , with 
Mr = ^^"^^ the WIMP-nucleon reduced mass, nir = "^^.'"^ the WIMP-nucleus reduced mass, 
the WIMP mass, mjsf the nucleus mass, and A the atomic weight. F is the nuclear form 
factor; in the following analysis the Woods-Saxon form factor |33] will be used (a more complete 
discussion can be found in reference |44j ) . Let us note that we are assuming identical WIMP-proton 
and WIMP-neutron cross-sections, and that we are ignoring the spin-dependent interactions. The 
integration over velocities is limited to those which can give place to a recoil energy Er, thus there 
is a minimal velocity given by fmin(-E'r) = 

A word of caution is also needed here to clarify the results we shall present. As it is obvious 
from equation (jlip . the sensitivity of a direct detection experiment depends strongly on the local 
DM density. The value po — 0.385 GeV cm~^ stated previously, refers to the overall density. If DM 
consists of multiple components, the differential event rate will depend on the partial density of the 
i-th. component (along with its velocity distribution). 

We already saw that the WMAP constraints can be satisfied in rather restricted regions of the 
parameter space. So, in most of the parameter space the model's relic density is larger than the 
ACDM one as inferred from WMAP, whereas in some cases it can also become smaller. In the case 
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ri^o h? > Odm h'^, we can speak of a 'hard exclusion' by experimental data, since the predicted relic 
density cannot be larger than the measured one. But in the case Q^o h? < f^DM things are more 
complicated. If neutralinos are not enough to explain the total DM relic density, nothing excludes 
them being only one of the components of the total DM density. In this case, the local density of 
neutralinos should be renormalised so as to account for this feature. In this sense, when computing 
whether a parameter space point is detectable or not, we should use the correct local density value 
(one could assume, for example, that the local density fraction scales as the relic density one). 

In the following, we shall be ignoring this point. We shall be computing sensitivity lines con- 
sidering the local density as being constant over the parameter space. In this spirit, the sensitivity 
lines should be read with a little caution. They can be read safely with respect to the regions where 
both WMAP bounds are fulfilled (the upper and the lower), whereas the reader should keep the 
previous remarks in mind for the bulk of the parameter space. When we say that a region of the 
parameter space is detectable, this corresponds literally to the WMAP fully compatible regions. 

4.2 A XENON-like experiment 

The XENON experiment aims at the direct detection of dark matter via its elastic scattering off 
xenon nuclei. It allows the simultaneous measurement of direct scintillation in the liquid and of 
ionization, via proportional scintillation in the gas. In this way, XENON discriminates signal from 
background for a nuclear recoil energy as small as 4.5 keV. Currently, the collaboration is working 
with a 170 kg detector, but the final project is a machine containing 1 ton of xenon. 

In our study, following reference |45| we will always consider 7 energy bins between 4 and 30 keV. 
We could take into account non-zero background using simulations of the recoil spectra of neutrons 
in our analysis, and this would significantly degrade the sensitivity of the detector. However, this 
would involve a much more detailed study of the detector components (shielding, etc.), and we will 
not carry it out. In that sense, our results will be the most optimistic ones. Comprehensive studies 
about the influence of astrophysical and background assumptions can be found in references [461 ITT] . 
Furthermore, we examine three 'benchmark' experimental setups, assuming exposures e = 30, 300 
and 3000 kg-year, which could correspond e.g. to a detector with 1 ton of xenon and 11 days, 4 
months or 3 years of data acquisition, respectively. 

Let us note that other promising direct dark matter detection experiments such as SuperCDMS 
|48) and LUX |49| . should give rise to sensitivities comparable to XENON's one. 

4.3 Results 

4.3.1 Correlated stop-slepton masses 

Figure [1] shows the exclusion lines (black lines) for exposures e = 30, 300 and 3000 kg-year, on 
the [ttt-o, Tni/2] parameter space, for all other parameters as defined in section [2T1 The first-row 
plots correspond to plain mSUGRA scenarios whereas the second and third to the 'mSUGRA- 
like' benchmark, with the ei and (ei,e2) parameters turned on respectively. The plots on the left 
correspond to a choice tan /3 = 3 whereas the right-hand side ones to tan /3 = 10. These curves 
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reflect the XENON sensitivity and represent its ability to test and exclude diflFerent regions of the 
mSUGRA and BMSSM relevant model at 95% CL: all points lying below the lines are detectable. 
We note that when some line is absent, this means that the whole parameter space can be probed 
for the corresponding exposure. Furthermore, the regions in orange (light gray) or blue (dark gray) 
are excluded due to the fact that the LSP is the stau or because of the null searches for charginos 
and sleptons at LEP. For large tan/3, an important fraction of the [mo, ?7ii/2] plane, corresponding 
to the region above the violet line, generates an unstable vacuum and is then excluded. Let us note 
that the introduction of £2 alleviates the vacuum stability constraint [11], and slightly increases the 
Higgs mass. 

As a general rule, the detection prospects are maximised for low values of the ttiq and 
parameters. For higher mo values, the masses of the squarks in the internal propagators increase, 
penalising the scattering cross-section. In the same way, the increase of mi/2 augments the WIMP 
mass and leads to a deterioration of the detection perspectives. On the other hand, the region of 
low mx/2 is also preferred because in that case the lightest neutralino is a mixed bino-higgsino state, 
favouring the Xi ~ Xi ~ ^ and Xi ~ Xi ~ H couplings, and therefore the scattering cross-section. 
Let us recall that a pure higgsino or a pure gaugino state does not couple to the Higgs bosons. 
On the other hand, the detection prospects are also maximised for low values of tan /3. For large 
values, besides the increase of the lightest Higgs boson mass, the coupling of the latter to a Xi pair 
decreases significantly because it is proportional to sin 2/3, for [//[ ^ Mi. 

The introduction of the NR operators gives rise to an important deterioration of the detection 
prospects. The main effect enters via the important increase in the lightest CP-even Higgs mass. 
This behavior is attenuated for larger values of tan/3; the corrections to the Higgs masses being 
suppressed by 1/ tan /3 (see e.g. reference 0). Moreover, as in the bulk region the lightest neutralino 
is mostly a bino-like state, the impact on its couplings with Higgs bosons is marginal. Nevertheless, 
let us emphasize again that this deterioration is relative, since we are comparing with a plain MSSM, 
which is already excluded because of the light Higgs mass. 

Concerning the plots in figure [H a further remark that can be made is that, even for low expo- 
sures, a sizable amount of the parameter space can be probed. The experiment will be particularly 
sensible to low values of mi/2- However, larger exposures could be able to explore almost the whole 
parameter space taken into account. Let us emphasize that, in general terms, the best detection 
prospects correspond to low values for tan/3 and for the couplings ej. 

4.3.2 Light stops, heavy sleptons 

Figure [2] shows the exclusion lines for XENON with exposures e = 30, 300 and 3000 kg-year, on 
the [Ml, //] parameter space, with the other parameters as defined in section fI7I\ for tan/3 = 3 (left 
panel) and 10 (right panel). Here again, the experiment will be sensitive to the regions below the 
contours. It can be seen that, in general, the detection prospects are maximised for low values for 
the Ml and/or the /i parameters, corresponding to a light Xv On the other hand, the scattering 
cross-section is enhanced near the region Mi ~ /i, where the lightest neutralino is a mixed bino- 
higgsino state, favouring the Xi~Xi~ ^ and Xi~Xi~ H couplings. Again, the detection prospects 
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Figure 1: Regions in the [niQ, mi/2] plane that can be detected by XENON using exposures e = 30, 
300 and 3000 kg-year, for our mSUGRA-hke scenario. The black lines depict the detectability 
regions: the area below the lines can be probed. Whenever a line is absent, this means that the 
whole parameter space can be tested by the experiment. The blue and orange regions depict the 
areas that are excluded by direct LEP chargino searches and the requirement for a neutralino LSP 
respectively. The area above the violet line is excluded by the metastable vacuum constraint. 
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Figure 2: Regions in the [Mi, fj] plane that can be detected by XENON for the scenario with hght 
stops and heavy sleptons. The black lines depict the detectability regions for the corresponding 
XENON detector with exposures e = 30, 300 and 3000 kg-year: the area below the lines can be 
probed. Whenever a line is absent, this means that the whole parameter space can be tested by the 
experiment. The blue and orange regions depict the areas that are excluded by direct LEP chargino 
searches and the requirement for a neutralino LSP respectively. The areas above the violet lines are 
excluded by the metastable vacuum constraints. 
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are also improved for low values of tan /3, mostly because the coupling between the LSP and the 
Higgs bosons is suppressed by the factor sin 2/3, for |^| ^ Mi. Let us note that the first line of the 
figure (corresponding to the case without the NR operators), besides being excluded by the Higgs 
mass, is partially ruled out by XENONIO 05] and CDMS \5Ql searches. 

When introducing the dimension 5 operators the detection prospects deteriorate, in a similar 
way as in the last subsection, because of the rise of the lightest Higgs mass. Furthermore, because 
of a suppression of the Xi ~ Xi ~ h coupling. The latter effect is very accentuated in the region 
where the LSP is liiggsino-like. It is interesting however to notice that, for the case of large tan/3, 
almost the whole area that could not be tested for e = 30 kg-year is already excluded by the 
vacuum stability constraint (i.e. the region above the violet line). Moreover, in this case the whole 
parameter space evades the current constraints from direct detection. 

Let us emphasise that this scenario offers moreover exceptionally good detection perspectives. 
Even with middle exposures, XENON will be able to detect dark matter in the whole region for all 
three benchmarks and in the two models. 



5 Gamma-rays from the Galactic Center 
5.1 Differential event rate 

The differential flux of gamma-rays generated from dark matter annihilations and coming from a 
direction forming an angle ^ with respect to the galactic center (GC) is 

°^^Xi "'^y -'line of sight 

where the discrete sum is over all dark matter annihilation channels, dN^/dE^ is the differential 
gamma-ray yield of SM particles into photons, {av) is the total self-annihilation cross-section 
averaged over its velocity distribution, p is the dark matter density profile, r is the distance from 
the GC and Bri is the branching ratio of annihilation into the i-th final state. The decay of SM 
particles into gammas has been calculated with PYTHIA [51]. The integration is performed along 
the line of sight from the observation point towards the GC. 

It is customary to rewrite equation (jl2p introducing the dimensionless quantity J, which depends 
only on the dark matter distribution (for an explicit calculation see, e.g. [521 [53] ): 

•^^ = 77^ /. p{r{l,^l^)f dl, (13) 

^0 Po J line of sight 

where Rq ~ 8.5 kpc is the distance from the Sun to the GC. After having averaged over a solid 
angle, Afi, the gamma-ray flux can now be expressed as 

<^y{Ey) = 1.55 • 10~^3cm~2s"^GeV~^sr-^ 

■ E^(,^)(^)^^(-)-. 
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The value of J{AQ)AQ depends crucially on the dark matter distribution. There are various ways 
to parametrize the density profile \54\ [5^1 156] . The most usual one is 

po[l + War] (Z^-^)/" 
' {r/RoV [1 + (r/a)°](/5-^)/" ' ^ ' 

where a is a characteristic length. 

There has been quite some controversy on the values for the (a, /3, 7) parameters. Some N-body 
simulations suggested highly cusped inner regions for the galactic halo |57) , whereas others predicted 
more moderate 7 values (the basic parameter determining the inner slope of the profile) |58j. The 
recent Via Lactea II simulation |59| seems to partly verify earlier results by Navarro et al. and finds 
their results as being well reproduced by (a, /3, 7) = (1, 3, 1). On the other hand, the Aquarius 
project simulation results |60) seem to be better reproduced by a different parametrization, the 
so-called Einasto profile |61) : 



pir) = exp 



2 // r . 

1 

a \ \r 



a = 0.17, (16) 



which does not demonstrate this effect of cuspyness in the inner galactic region. Here = 20 kpc 
is a characteristic length, while ps is a normalization factor, which we fix so as to reproduce the 
solar dark matter density. This choice yields ps ~ 0.0783 GeV cm^'^. 

In most numerical simulations, among which the Via Lactea II and Aquarius, it is a common 
simplification that the effect of baryons is not taken into account. It has however been pointed out 
that in the presence of baryons there can be adiabatic collapse phenomena taking place near the 
galactic center (see, e.g. references |62|I63)). something which could severely influence the innermost 
regions of the DM halo, leading to profiles much more cusped than the ones usually predicted, and 
hence enhancing the relevant fluxes by important factors. We thus also consider a profile which, 
starting from the NFW one, tries to take into account such effects, leading to an enhanced inner 
slope; this profile with adiabatic compression is denoted NFWc. The parameters relevant for the 
models under discussion can be seen in table [1] 

It is worth noticing here that we are neglecting the effect of clumpyness, even though other 
studies showed that, depending upon assumptions on the clumps' distribution, in principle an 
enhancement by a factor 2 to 10 is possible [641159) . In this respect, the following predictions on the 
gamma-ray flux from the galactic center are conservative. This shall also be the case in the following 
analysis regarding the detection capacity of the model in the positron and antiproton channels. 

As a flnal remark, let us repeat the word of caution already mentioned for direct detection. In 
order to draw sensitivity lines, we shall be considering a one-particle DM with the aforementioned 
density proflles irrespectively of the relic density inferred from the model. Hence, once again, these 
lines should be read with respect to the WMAP-compatible regions. 

5.2 Modeling the galactic center background 



HESS |65| has measured the gamma-ray spectrum of a very bright point-like source very close to 
the galactic center in the range of energy [160 GeV-10 TeV]. The collaboration claims that the 
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a [kpc] 


a 


/3 


7 


J(3 • 10"^ sr) 


Einasto 










6.07 • 10^ 


NFW 


20 


1.0 


3.0 


1.0 


8.29 • 10^ 


NFWc 


20 


0.8 


2.7 


1.45 


5.73 • 10^ 



Table 1: Einasto, NFW and NFWc density profiles with the corresponding parameters, and values of J(Af2). The 
latter has been computed by means of a VEGAS Monte-Carlo integration algorithm, imposing a constant density for 
r < 10~^kpc so as to avoid divergences appearing in the NFW-like profiles. 

data are fitted by a power-law 

<^™f (^)=i^0i?^eV' (17) 

with a spectral index a = 2.21 ± 0.09 and Fq = (2.50 ± 0.21) • 10"® m'^ s"^ TeY~\ The data were 
taken during the second phase of measurements (July-August, 2003) with a of 0.6 per degree 
of freedom. Because of the constant slope power-law observed by HESS, it turns out possible 
but difficult to conciliate such a spectrum with a signal from dark matter annihilation \63\ 166) . 
Indeed, final particles (quarks, leptons or gauge bosons) produced through annihilations give rise to 
a spectrum with a continuously changing slope. Several astrophysical models have been proposed in 
order to match the HESS data [67] , In the present study we consider the astrophysical background 
for gamma-ray detection as the one extrapolated from the HESS data with a continuous power-law 
over the energy range of interest (~ 1 - 300 GeV). 

EGRET |68) reported the presence of a bright gamma-ray source at energies below 10 GeV, 
which exceeds by far the HESS aforementioned extrapolation. However, this source seems not to 
be confirmed by the recent Fermi collaboration data |69| . We shall hence not take into account this 
point source. 

Finally, we will consider the diffuse background of gamma rays in the region surrounding the 
galactic center. We will describe the spectrum of the background using the HESS observation from 
the Galactic Center Ridge [70], which can be described by 

(PhkgiE) = 1.1 • 10"^ E^'^^^ GeV^^cm-^s-isr-i . (18) 

In our analysis, we shall consider a solid angle of observation around the galactic center (AQ = 
3 • 10~^ sr) and the energy region between 1 and 300 GeV. 

5.3 The Fermi experiment 

The space-based gamma-ray telescope Fermi |7H I72| was launched in June 2008 for a five-year 
mission. It performs an all-sky survey covering a large energy range (~ 30 MeV - 300 GeV). 
With an effective area and angular resolution on the order of 10^ cm^ and 0.1° (Afi ~ 3 • 10~^ 
sr) respectively, Fermi will be able to point and analyze the inner center of the Milky Way (~ 7 
pc). Concerning the requested condition on the for a signal discovery, we have used an analysis 
similar to the one considered in the case of direct detection as defined in section [3l Additionally, 
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we take into account a five-year mission run, and an energy range extending up to 300 GeV, with 
20 logarithmically evenly spaced bins. 

5.4 Results 

5.4.1 Correlated stop-slepton masses 

In figure [3] we present the detectability regions for a 5- year run of the Fermi experiment and for 
3 different halo profiles presented in the literature, Einasto, NFW and NFWc, in the [rriQ, mi/2] 
parameter space. Fermi will be sensitive to the regions below the contours and, for tan/3 = 3, to 
the area inside the blob. The detection prospects are maximised for low values of the niQ and 
parameters. For higher thq values, the masses of the squarks increase, penalising the annihilation 
cross-section. However, the growth of mi/2 gives rise to the opening of some relevant production 
channels, after passing some thresholds, increasing significantly the {(Jv). The first one corresponds 
to a light neutralino, with mass ~ mzf^ ("^1/2 ~ 130 GeV). In that case the annihilation is 
done via the s-channel exchange of a real Z boson, decaying in hadrons (~ 70%), neutrinos (~ 20%) 
and charged leptons (~ 10%). The second threshold appears for rriy. ~ mw ~ 220 GeV). The 

annihilation cross-section is enhanced by the opening of the production channel of two real 
bosons in the final state. This process takes place solely through chargino exchange, since both Z 
and Higgs bosons exchange are suppressed by taking the limit f — )• 0. The last threshold corresponds 
to the opening of the channel XiXi ~^ ("I1/2 ~ 400 GeV). The diagrams involved in such a process 
contain contributions from t- and n-channel exchange of stops, and from s-channel exchange of Z's 
and pseudoscalar Higgs bosons. The aforementioned threshold appears as a particular feature on 
the left-hand side plots: An isolated detectable region for mi/2 ~ 400-500 GeV and tjiq < 300 GeV 
corresponding to the annihilation into a pair of real top quarks. 

Larger values for the annihilation cross-section can be reached for higher values of tan /3. In 
that case, the production process of a pair of down-type quarks (in particular bb pairs) and charged 
leptons, dominates the total cross-section. In fact, the diagrams containing exchanges of a pseu- 
doscalar Higgs boson or a sfermion are enhanced by factors tan (3 and 1 / cos (3 respectively. On the 
other hand, for high values of tan /3, the channels corresponding to the annihilations into W~^W~ 
and vanish. The first because of the reduction of the coupling Xi ~ xt ~ the second because 
of important destructive interference between diagrams containing the exchange of a Z boson and 
stops. 

For the present scenario, the introduction of the NR operators gives rise to a very mild signature. 
Actually, as in almost the whole parameter space the lightest neutralino is bino-like, its couplings 
do not vary drastically. Moreover, the increment in the Higgs masses has a small impact on the 
(av) factor. For indirect detection prospects, the main effect corresponds to a slight increase in the 
LSP mass. Let us emphasize on the fact that, however, the detectable regions are in the BMSSM 
case more cosmologically relevant than in the corresponding plain MSSM one. 

Concerning figure [3l let us note that the only astrophysical setup in which some useful infor- 
mation can be extracted is the NFWc one. This means that in this scenario, in order to have 
some positive detection in the 7-ray channel, there should exist some important enhancement of 
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Figure 3: Regions in the [mo, 'm'i/2] plane that can be detected by the Fermi satehite mission for 
our mSUGRA-Hke scenario. The black lines depict the detectability regions for the corresponding 
halo profile assumptions and 5 years of data acquisition: the area below and on the left of the lines 
can be probed. The same applies to the top-resonance blob at mi/2 ~ 450 GeV appearing on the 
left-hand side plots. For NFW and Einasto profiles, the model could not be tested. 
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the signal by some astrophysical mechanism (as the adiabatic contraction mechanism invoked in 
this case). We note that, and this wih be different from the case of antimatter signals, there is 
however no important constraint on astrophysical boosts from the Galactic Center. Gamma-ray 
detection does not rely, as is the case for positrons that we shall examine in section [6l that much on 
local phenomena. In this respect, the NFWc results can be characterized as optimistic (it has been 
pointed out that even by changing the gravitational collapse conditions, the results can get even 
more pessimistic in the case, e.g., of a binary black hole formation in the GC), but not excluded. 

5.4.2 Light stops, heavy sleptons 

Figure 2] presents the results for the second scenario with light stops and heavy sleptons. The 
experiment will be sensitive to the regions below/on the right of the contours. Again, the detection 
prospects are maximised for low values of the Mi and parameters, corresponding to light WIMPs. 
However, the growth of any of the latter parameters gives rise to the opening of some production 
channels, enhancing significantly the (crv). The first one appears for ~ 'mz/'2 and corresponds 
to the s-channel exchange of a real Z boson. The second one concerns the production channel of 
two real W bosons. Let us note that in this scenario the neutralino LSP can be as heavy as ~ 110 
GeV, implying that the annihilation into a pair of top quarks is never kinematically allowed. On the 
other hand, the region where Mi ^ n is highly favored for indirect detection due to the fact that 
the LSP is higgsino-like, maximising its coupling to the Z boson. Let us recall that the Z boson 
does not couple to a pure gaugino-like neutralino. 

Large values for the annihilation cross-section can be reached for high values of tan /3, mainly 
because of the enhanced production of bb pairs. On the other hand, for high values of tan (3, the 
threshold corresponding to the opening of the annihilation into W~^W~ is suppressed or enhanced 
for /i ^ Ml or /i ^ Mi respectively, due to the dependence of the Xi ~ xf ~ coupling on the 
texture of the LSP. 

For the present scenario, the introduction of the NR operators gives rise to an important increase 
of the Xi — Xi ~ ^ coupling when fi > Mi , and therefore to a boost in the annihilation into fermion 
pairs. On the other hand, as the Higgs boson h becomes heavier, the processes giving rise to the 
final state hZ get kinematically closed. 

In the case presented in figure there is a positive detection for all three halo profiles; however, 
the regions that can be probed for either the NFW or the Einasto cases are cosmologically irrelevant. 

In fact, they could give rise to a positive detection near the Z-funnel and in the region where the 
LSP is a higgsino state {Mi > 150 GeV); nevertheless the first is already excluded by LEP (at least 
for minimal scenarios) and the second generates too small a dark matter relic density, below the 
WMAP limits. On the other hand, the profile NFWc could test a large amount of the parameter 
space we examine, particularly for high values of tan (3. Only the Higgs peak and the regions with 
a heavy LSP escape from detection. 
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Figure 4: Regions in the [Mi, fj] plane that can be detected by the Fermi satehite mission for our 
scenario with hght stops and heavy sleptons. The black lines depict the detectability regions for the 
corresponding halo profile assumptions and 5 years of data acquisition: the area below the lines can 
be probed. The same applies to lines forming closed regions with respect to the axes, as is the case 
for the NFW and NFWc profiles: the parameter space points lying in the interior of these regions 
yield signals that are detectable. 
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6 Positrons 



Positrons, as antiprotons that we shall describe in the following, present the complication of pro- 
pagating throughout the Intergalactic Medium. Numerous treatments, making a different set of 
assumptions and simplifications have been presented in the literature |73| [TH [75] and applied in the 
case of the MSSM |76| \77\ \78\ I79| . Here we shall use the two-zone diffusion model and its semi- 
analytical solution as described in reference [75]. In this model, positron and antiproton propagation 
takes place in a cylindrical region (Diffusive Zone, DZ) around the galactic center of half thickness 
L; the propagating particles being free to escape the region, a case in which they are simply lost. 
The propagation is described by a diffusion-convection-reacceleration equation: 

^t^P + ^,iVc^P)-ViKVi,)-^E[biE)^P + KEEiE)^E^]=q, (19) 

where ip = dn/dE is the space-energy density of the positrons or antiprotons, h{E) is the energy 
loss rate, q is the source term, ~ (5 — 15) km/s is the convective wind velocity wiping away the 
positrons or antiprotons from the galactic plane, 

K{E) = K,^I^^^ (20) 

is the diffusion coefficient, with (3 being the particle's velocity, Kq the diffusion constant, a a 
constant slope, E the kinetic energy (for positrons in practice the total one), Eq a reference energy 
(which we take to be 1 GeV), and 

= (21) 
is a coefficient describing reacceleration processes. 

This is the master equation governing the propagation of cosmic rays throughout the galactic 
medium, which we shall be employing in the following for both positrons and antiprotons, with 
the details varying, of course, according to each species (e"*" or p). 

6.1 Differential event rate 

It has been pointed out in the literature that the convective wind as well as reacceleration processes 
can be quite safely neglected for the case of positrons [80 j (at least to a level sufficient for our 
purposes). On the other hand, energy losses should be taken into account, the most important 
contributions coming from synchrotron radiation or inverse Compton scattering on stellar light and 
CMB photons. To account for these two processes, we write the energy loss rate as: 

b{E) = , (22) 

Eq te 

where E is the positron energy and te = lO^^s is the characteristic energy-loss time. We are then 
left with the following equation: 

dti^ - V [K{x, E) W] - Be [b{E) V^] = q{x, E) , (23) 
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L [kpc] 


[kpcVMyr] 


a 


MIN 


1 


0.00595 


0.55 


MED 


4 


0.0112 


0.70 


MAX 


15 


0.0765 


0.46 



Table 2: Values of propagation parameters widely used in the literature and that roughly provide minimal and 
mciximal positron fluxes, or constitute the best fit to the B/C data. 



where K is the space diffusion coefficient -steady state is assumed. This coefficient is taken to be 
constant in space but depends on the energy as 



E 



a 



Here the diffusion constant, Kq, and the spectral index, a, are propagation parameters. Regarding 
the propagation parameters L, Ko and a, we take their values from the commonly used MIN, MAX 
and MED models -see table [2l The former two models correspond to the minimal and maximal 
positron fluxes that are compatible with the B/C data [81j. The MED model, on the other hand, 
corresponds to the parameters that best fit the B/C data. 

The master equation for positron propagation (equation ()23p ) gets simplified to its final form 

^oe"VV + ^ +9 = 0, (25) 

where e = E/Eq. This is the expression we solve to calculate the effects of positron propagation on 
a signal produced at some point in the galaxy. 

The resulting positron fiux from DM annihilations can be written as (see reference [73l I75| for 
details) 

where the detection and the production energy are denoted respectively by E and Eg, Xq is the solar 
position with respect to the GC and /3g+ is the positron velocity. /{Eg) is the production spectrum 
for positrons, f{Es) = "^^dN^^/dEg, with i running over all possible annihilation channels. The 
diffusion length, A/j, is defined by 

'a—l ^a- 1 



Al, = 4A„r.(^-^j. (27) 

The so-called halo function, /, contains all the dependence on astropliysical factors. It is given by 

i{\D)= I d^XsG{xQ,E^Xs,Es) (4^) , (28) 
Jdz \p{xq)J 

where the integral is performed over the diffusion zone. The modified Green function G is in its 
turn defined by 

G = ^ g-(r0-r,)V(4i^o^) Y (29) 

AirKof 
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L2 



and f 



t - t. 



with V depending on the value of the characteristic parameter ( — 
te [(^"""'^/(l ~ Oi)) ~ (^""^/(l ~ ct))] ■ When C > 1 -when the diffusion time is small- boundary 
conditions can be ignored and the propagation equation can be treated as a ID Schrodinger equa- 
tion. In that case 



V 



1 



= exp 



When C, is small this approximation no longer holds but we can express V as 



where 



V 



oo ^ 



n=l 



sin[A:^(L 



n — 



vr 
I 



(30) 

(31) 

(32) 
(33) 



A„ = Ko k"^ and = Kq {k'^)"^ . 

The advantage of this method is that the halo function I{\d) can be calculated (and either 
tabulated or fitted) just once as a function of the diffusion length and then be easily used for 
performing parameter space scans which, as in our case, can be rather large. 

Let us however note, and this shall also be the case for antiprotons, that this method has 
been proved to have a limited validity. More specifically, at energies below 10 GeV, some of the 
assumptions and simplifications that we made no longer hold. For this reason we shall limit ourselves 
at energies > 10 GeV. 

6.2 The background 

In the conventional background model [82], positrons are produced in the interactions between 
cosmic-ray nuclei and the interstellar medium (ISM). This model is however not compatible with 
the recent data from PAMELA [Blj and Fermi-LAT [Si] . 

Even after taking into account several possible uncertainties due to cosmic ray propagation, the 
data reveals a clear excess over this background at high energies, S > 10 GeV. Hence, a new source 
of high energy positrons is necessary to explain the data. It also turns out to be exceedingly difficult 
to reconcile the observed excess with dark matter annihilations. A large number of models trying 
to do so have been proposed in the literature, most of which turn out to be in confiict with other 
observational data (see, e.g. the interesting treatments in references |85l l86l [871 EHl |89l |90l [91] [92| 
[931 IMl [951 [961 [97]). 

Although obviously dark matter annihilation could, in principle, account at least for some part 
of the observed excess, we feel that the most conservative choice is to consider the whole PAMELA 
signal as been due to some -yet unknown- astrophysical process, such as a pulsar [981 [99l 1100) and 
treat it as a background for the oncoming AMS-02 experiment |101) . We shall use as a background 
the absolute positron fiux that can be obtained through a combination of the PAMELA and Fermi 
data. 
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6.3 AMS-02 and positron detection 

The scheduled AMS-02 mission aims at the detection and measurement of cosmic-ray fluxes (and 
also 7-rays) coming from various sources. Among these sources could, of course, be dark matter 
annihilations. 

In the case of positrons, AMS-02 will be able to measure the spectrum of positrons with an 
average geometrical acceptance of 0.042 m^ sr in the energy range above 4 GeV |102) . In our study 
we consider a 3-year run, which is actually the collaboration's nominal run time, and an energy 
range extending up to 300 GeV, with 20 logarithmically evenly spaced bins. 

6.4 Results 

6.4.1 Correlated stop-slepton masses 

The results concerning the detectability perspectives for the mSUGRA-like scenario in the positron 
detection channel are quite pessimistic. In fact, since the PAMELA and Fermi measurements 
mentioned above, and according to our conservative treatment, the main issue in the positron 
channel is an extreme domination of all measurements by a large background severely obscuring the 
signal. 

Obviously, one could invoke large boost factors of an astrophysical nature as was the case in 
the first efforts to explain the PAMELA excess through dark matter annihilations, a case in which 
a larger portion of the parameter space would be visible. However, it has been pointed out that it 
is highly unlikely to expect large boost factors due, e.g., to substructures in the halo |103) . In this 
respect, if we assume a maximal clump-due signal enhancement by a factor ~ 10, the only hope 
for positive detection of a non-LEP excluded area might come for the bulk region, as it is the only 
one lying at the limits of detectability. For the sake of brevity, we omit the relevant plots for the 
mSUGRA-like benchmark, since no point of the parameter space can be tested. 

6.4.2 Light stops, heavy sleptons 

In figure [5] we present the detection perspectives in the positron channel for our scenario with light 
stops and heavy sleptons. The detectable parameter space regions lie within the zones delimited 
by the black lines for the three propagation models: the oval-shaped blobs as well as the banana- 
shaped ones. Once again, we notice the general features already present in the 7-ray channel. The 
regions giving rise to a positive detection lie within the zone where the LSP is a higgsino-like state, 
with mass > mw, in order to have the final state W~^W~ kinematically available. This region 
in general does not fulfill the WMAP limit. However, and this is a novel feature of the BMSSM, 
with both ei and €2 couplings turned on, a small region of the mixed higgsino-bino regime can be 
detected for the MAX (and even the MED) propagation model. As we pointed out before, in this 
regime the total annihilation cross-section can be quite significantly enhanced, leading to better 
detection perspectives. 



22 



tan p=3, E,=0, E2=0 



tan p=10, Ei=0, E2=0 




400 



> 
O 



> 

o 



500 



400 



300 



100 150 
M, [GeV] 

tan p=3, E,=-0.1, 



100 150 
M, [GeV] 

tan p=3, £,=-0.1, E2=0.05 



250 




250 




250 



> 

o 



> 
o 



> 
o 




WMAP 

r, Lsp 

LEP 

e"*" Max 
e^Med 
e"*" Min 



100 150 
M, [GeV] 

tan|3=10, £,=-0.1, 



250 




100 150 
M, [GeV] 

tanp=10, E,=-0.1, £,=0.05 




WMAP 
r, LSP 
LEP 
V. stability 
e Max 
e^Med 
e"*" Min 



250 



Figure 5: Regions in the [Mi, fj] plane that can be detected by a 3-year run of the AMS-02 satelhte 
mission for the scenario with hght stops and heavy sleptons, in the positron channeh The black lines 
depict the detectability regions for the 3 considered propagation models, MIN, MED and MAX: 
the parameter space points lying within the regions delimited by the black lines can be probed, 
assuming the corresponding propagation models. Part of the mixed bino-higgsino region, as well as 
(marginally) some part of the Z funnel region can be probed. 
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7 Antiprotons 

7.1 Differential Event Rate 



Equation ()19p also governs antiproton propagation in the galactic medium. Nevertheless, the do- 
minant processes vary significantly with respect to the positron case. More precisely, all energy 
redistributions in the initial (injection) spectrum -energy losses, reacceleration, as well as 'tertiary' 
contributions (i.e. contributions from secondary antiprotons produced upon inelastic scattering 
with the interstellar medium)- can be ignored. Whether these redistributions are important or not 
depends mainly on the antiproton energy. For GeV energies, the results may deviate up to 50% 
from those obtained with the (more complete) Bessel function .reatme^a But for energies around 
10 GeV, the accuracy of the method improves dramatically, yielding essentially indistinguishable 
results at slightly higher energies. Since the p energy region we shall consider begins at 10 GeV, we 
can safely use this simplified approach. 

Let us denote by F^"" = X^ism '^ism v cr^i^^ the destruction rate of antiprotons in the interstellar 
medium, where ISM = H and He, ti-ism is the average number density of ISM in the galactic disk, 
V is the antiproton velocity and o"p""g^ is the p — ISM annihilation cross-section. Implementing the 
aforementioned simplifications, the transport equation for a point source (which actually defines 
the propagator G) is: 

-i^v + y,|- + 2/ir|^^5(z 

oz ^ 



G = 6[r-r'] , (34) 



with h = 100 pc being the half-thickness of the galactic disc. The antiproton propagator at the 
solar position can then be written (in cylindrical coordinates) as 

G^{r, z) = 1^ Ko (rVfc2+A;2) sin(/c„ L) sin(A:„ (L - z)) , (35) 

n=0 

where Kq is a modified Bessel function of the second kind and 

sin(A;„L) cos(A;„L) 
Cn = 1 ' ^^^) 

K = VJ(2K), (37) 

kd = IhT^'^/K + 2k^ . (38) 

kn is obtained as the solution of the equation 

nir — kn L — arctan{2kn/k(i) = 0, riGN. (39) 



Then, in order to compute the fiux expected on earth, we should convolute the Green function 
with the source distribution q{x,E). For dark matter annihilations in the galactic halo, the source 
term is given by 



^In reference [104] a comparison between the two methods can be found (see figure 2). 
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L [kpc] 


Ko [kpcVMyr] 


a 


Vc [km/s] 


MIN 


1 


0.0016 


0.85 


13.5 


MED 


4 


0.0112 


0.70 


12.0 


MAX 


15 


0.0765 


0.46 


5.0 



Table 3: Values of propagation parameters widely used in the literature and that provide minimal and maximal 
antiproton fluxes, or constitute the best fit to the B/C data. 

where the index i runs over all possible annihilation final states. As in the previous cases, the decay 
of SM particles into antiprotons has been calculated with PYTHIA |51) . Regarding the distribution of 
dark matter in the Galaxy, p{x), we assume a NFW profile. The final expression for the antiproton 
flux on the Earth takes the form 

where none of the integrated quantities depends on the antiproton energy. Once again, the integral 
in equation ()4ip . which we compute using a VEGAS Monte-Carlo algorithm, needs to be calculated 
only once for each value of the injection energy, which is actually the same as the detection energy. 

Regarding the propagation parameters L, Kq, a, and V^, we take their values from the well- 
established MIN, MAX and MED models -see table |3j The former two models correspond to the 
minimal and maximal antiprotons fluxes that are compatible with the B/C data. The MED model, 
on the other hand, corresponds to the parameters that best fit the B/C data. 



7.2 The background 

Contrary to the positron case, the most well-known treatment of the astrophysical antiproton back- 
ground by Strong and Moskalenko [105] . seems to be compatible with the antiproton data |106) 
from the PAMELA experiment. 

In order to parametrize our background, we borrow the simple fit performed by Cirelli et al. 
|107| . which provides a sufficiently good fit for our purposes, to theoretical predictions but also to the 
recent PAMELA data. We pay special attention at reproducing the good background normalization 
at low energies, so as to stay as close as possible to the PAMELA measurements. 

7.3 AMS-02 and antiproton detection 

In the case of antiprotons, AMS-02 will be able to measure [102] the corresponding fluxes with an 
average geometrical acceptance of 330 cm^ sr above 16 GeV and up to 300 GeV. In our study we 
again consider a 3-year run and the mentioned energy range divided into 20 logarithmically evenly 
spaced energy bins. 
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7.4 Results 

7.4.1 Correlated stop-slepton masses 

In figure [6] we present our results for the detectability of the BMSSM in comparison to the MSSM by 
the AMS-02 experiment for the antiproton channel. The detectable regions lie below the black lines. 
In the case tan (3 = 3, the experiment is not sensitive to any point in the parameter space satisfying 

tan|3=10, e,=0, £2=0 




100 300 500 700 900 
mo [GeV] 

tanp=10, £[=-0.1, £2=0 tanp=10, E,=-0.1, £2=0.05 




mo [GeV] mo [GeV] 

Figure 6: Regions in the [mo, ?tT'1/2] plane that can be detected by a 3-year run of the AMS-02 
satellite mission for our mSUGRA-like scenario in the antiproton channel. The black lines depict 
the detectability regions for the 3 considered propagation models: the area delimited by the axes 
and the black lines can be probed for the corresponding propagation model (i.e. the region towards 
the lower left corner in each plot). 

also the collider constraints (and, hence, the corresponding results are once again omitted). A first 
remark here should concern the fact that the perspectives for antiproton detection are significantly 
ameliorated with respect to the corresponding positron ones, at least for large values of tan /3. 
This could, in some sense, seem quite strange, since the average antiproton yield coming from a DM 
annihilation is, in general, smaller than the positron one. However, we should take into account that 
these results reflect not only the behaviour of the models themselves, but also the interplay between 
the signal and the background. And, as it turns out the p channel is a signiflcantly low-background 
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one. 

Important areas of the viable parameter space are at the hmits of detectability: the bulk region, 
but also, for some cases, the Higgs funnel. Now, as we stressed out before, the possible enhancements 
due for example to substructures are quite constrained. Given however that some regions are 
marginally out of reach, it would not be impossible to state that even small boosts could render 
important (in a qualitative sense, due to their cosmological relevance) regions of the parameter 
space detectable by AMS-02. 

7.4.2 Light stops, heavy sleptons 

Figure [7] presents the results for antiprotons and for the second scenario under consideration. AMS- 
02 will be able to probe the regions lying within the oval-like blobs and the banana-shaped regions 
delimited by the black contours and the axes. Once again, the BMSSM turns out to be more 
favorable for DM detection than the ordinary case of light stops and heavy sleptons without NR 
operators. Detectable regions fall either into the case of the higgsino-like neutralino regime, or in 
the low-mass Z funnel region. We point out that an important part of the area where the dark 
matter relic density is fulfilled via coannihilation with the lightest stop could also be tested. 

8 Conclusions 

In this paper, we presented detection perspectives for neutralino dark matter in the framework of 
the so-called BMSSM, where the superpotential is enriched by the addition of dimension-5 non- 
renormalizable operators in the Higgs sector. The main motivation to add non-renormalizable 
operators to the MSSM Higgs sector is to reduce the fine-tuning that is required by the lower bound 
on the Higgs mass. These new terms re-open the bulk region, already excluded because it gives rise 
to a too light Higgs boson. Focusing on the new available regions, we have studied four of the most 
popular detection modes: direct detection, as well as three different channels of indirect detection, 
7-rays from the galactic center, positrons and antiprotons. We placed ourselves in the framework 
of the experiments XENON, Fermi and AMS-02. 

According to our results, the most favourable detection mode seems to be, by far, the direct 
one. XENON, even for low exposures, can provide a sizable parameter space coverage and, in 
case of a negative detection, can exclude a significant part of the model. Moreover, regions of the 
supersymmetric parameter space that were favored by the WMAP constraint, but were excluded 
within the MSSM, become viable within the BMSSM and could be potentially detected. 

From the indirect detection point of view, we have shown that gamma-rays are the most favoured 
probe in the benchmarks of the BMSSM we examined. We stress nevertheless that this is the case 
assuming quite optimistic astrophysical conditions, namely an important collapse of dark matter 
in the galactic center due to the presence of baryons, yielding a significantly spiked inner region 
profile. We should also note again that we considered a 5-year Fermi mission, whereas a period of 
3 years was considered for the AMS-02 experiment. On the other hand, antiprotons can provide 
significant information without having to assume extreme conditions for the galactic medium: we 



27 



tan p=3. 



=0, Ei=0 



tan p=10, Ei=0, E2=0 



> 
O 



500 



400 



300 




400 



> 

o 



> 
o 



500 



400 



300 



100 150 
M, [GeV] 

tan p=3, E,=-0.1, 



250 




100 150 
M, [GeVJ 

tan p=3, £,=-0.1, £,=0.05 



250 




250 



> 

o 



> 
o 



> 
o 




100 150 
M, [GeV] 

tan|3=10, £,=-0.1, 



100 150 
M, [GeVJ 

tanp=10, £,=-0.1, £,=0.05 



WMAP 
t, LSP 
LEP 
p Max 
p Med 
p Min 



250 




WMAP 
t, LSP 
LEP 
V. stability 
p Max 
p Med 
p Min 



250 




WMAP 
r, LSP 
LEP 
V. stability 
p Max 
p Med 
p Min 



250 



Figure 7: Regions in the [A/i, fi] plane that can be detected by a 3-year run of the AMS-02 satelhte 
mission for the scenario with hght stops and heavy sleptons, in the antiproton channel. The black 
lines depict the detectability regions for the 3 considered propagation models: the areas delimited 
by the axes and the black lines can be probed for the corresponding propagation model (i.e. the 
regions towards the lower right side in each plot), as well as the areas delimited by closed lines. 



only used the standard propagation models used in numerous other analyses in the literature. 
However, the detection prospects could be enhanced taking into account DM substructures, which 
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could provide a -small- amelioration of the signals by at most one order of magnitude. For positrons, 
our treatment for the background seems to lead us to conclude that the PAMELA data, providing 
a large supplementary background with respect to previous estimates, seem to render a positive 
detection in this channel more difficult. 

As a final remark, it is interesting to note again that for the low tan (3 values we considered here, 
as we wanted to stick to the original motivation for the model (uplifting the Higgs mass without 
over-constraining the stop sector), direct detection seems dominant. It is quite well-known that 
increasing tan /3 usually tends to ameliorate the indirect detection prospects, contrary to direct 
detection ones. This is an element which can give us an idea of the complementarity of the various 
detection modes. As underlined in I91|. different detection techniques offer the possibility for 
a more complete parameter space coverage. Moreover, they can serve as a means for independent 
confirmation of discoveries and/or comparison with other constraints. Given the controversy that 
has been generated during the past few years on the nature of the various excesses that have 
occasionally been observed, we feel that combining information from various sources is an essential 
element in our searches. 
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